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Abstract 

A formalism based on a chiral quark model (xQM) approach complemented with a one-gluon 
exchange model, to take into account the breakdown of the SU(6) £3 0(3) symmetry, is presented. 
The configuration mixing of wave functions for nucleon and resonances are derived. With few 
adjustable parameters, differential cross-section and polarized beam asymmetry for the 7 p — > n p 
process are calculated and successfully compared with the data in the centre-of-mass energy range 



> • from threshold up to 2 GeV. The known resonances 5n(1535), 5n(1650), Pi 3 (1720), Di 3 (1520), 

and ^15(1680), as well as two new Su and D15 resonances are found to be dominant in the reaction 
mechanism. Besides, connections among the scattering amplitudes of the xQM approach and the 
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helicity amplitudes, as well as decay widths of resonances are established. Possible contributions 

00 ' 

from the so-called "missing resonances" are investigated and found to be negligible. 
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I. INTRODUCTION 



Electromagnetic production of mesons on the nucleon offers a great opportunity to deepen 
our understanding of the baryon resonances properties. In recent years, intensive experi- 
mental efforts have been devoted to the measurement of observables for the processes of 
pseudoscalar and vector mesons production, using electron and/or photon beam facilities. 

In the present work we investigate the reaction 7 p —>■ 77 p, in the range of centre-of- 
mass total energy from threshold up to W ~ 2 GeV, in order to interpret a large amount 
of high quality data released from various facilities, namely, differential cross-section data 
by the following collaborations: MAMI fl, CLAS Q, CB-ELSA Q, LNS-GeV-7 Q and 
GRAAL 5], polarized beam asymmetries by CB-ELS A/TAPS and GRAAL 

The copious set of data has motivated extensive theoretical investigations. Most of the 
available models are based on meson-nucleon degrees of freedom, in which the Feynman 
diagrammatic techniques are used, so that the transition amplitudes are Lorentz invariant. 
In recent years various advanced approaches have been developed, namely, the unitary isobar 



model of MAID [7], Geissen [8] and Bonn-Gatchina groups [9i| coupled-channel approaches, as 
well as the partial wave analysis of SAID [lo| . Those approaches have no explicit connection 
with QCD, and the number of parameters in the models increases with the number of 
resonances included in the models. 

Formalisms embodying the subnucleonic degrees of freedom are also being developed. 
Such a program has its genesis in the early works by Copley, Karl and Obryk H[ and Feyn- 
man, Kisslinger and Ravndal 12| in the pion photoproduction, who provided the first clear 
evidence of the underlying SU(6) ® 0(3) structure of the baryon spectrum. The subsequent 



works 
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141 ] in the framework of the constituent quark models concentrated mainly on 



the transition amplitudes and the baryon mass spectrum, predicting still undiscovered or 
"missing", resonances. However, those approaches did not investigate reaction mechanisms. 

In Ref. [15|] a comprehensive and unified approach to the pseudoscalar mesons photo- 
production, based on the low energy QCD Lagrangian [if]], is developed with the explicit 
quark degrees of freedom. This approach reduces drastically the number of free parameters, 
for example, within the exact SU(6) ® 0(3) symmetry, the reaction under investigation 
has only one free parameter, namely, r]NN coupling constant. However, that symmetry is 
broken and in order to take into account that effect, one free parameter per resonance was 
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introduced in previous calculations [l?j, Q]- Given that the configuration mixing among 
the 3-constituent quarks bound states is a consequence of the SU(6) (g> 0(3) symmetry 
breakdown, in the present work we use the one-gluon-exchange mechanism to generate the 
configuration mixing of the wave functions. In this approach, the number of parameters de- 
creases significantly. After the parameters are determined by fitting the data, we then study 
the contributions from the missing resonances (see e.g. Refs. IjJ, [2(J, [2l| ) . Besides, we give 



relations connecting the scattering amplitudes in our xQM approach to the photoexcitation 
helicity amplitudes and partial decay widths of resonances. Our approach offers also the 



opportunity of investigating new nucleon resonances^ 



reported in the literature 
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22 



23 



24 



26 



27 



or w 
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rich str ong indications have been 



29 



30 



3l|. 



The paper is organized as follows. In Section [TTl the theoretical content of our work is 
presented. Starting from a chiral effective Lagrangian, the CGLN amplitudes for the process 
7 p — > T] p are given within the SU(6) <S> 0(3) symmetry. Then the consequences of the 
breaking of that symmetry via configuration mixing in One-Gluon-Exchange (OGE) model 
is reported and helicity amplitudes of photon transition and meson decay partial widths 
of resonances are presented. The fitting procedure and numerical results for differential 
cross-section, polarized beam asymmetry, helicity amplitudes, and N* — > r] N partial 
decay width are reported and discussed in Section IHH where possible roles played by new / 
missing resonances are examined. Summary and conclusions are given in Section IIVI 



II. THEORETICAL FRAME 



In this Section we recall the content of a chiral constituent quark approach and relate 
it to the configuration mixing of constituent quarks states via a OGE model, generated 
by the breakdown of the SU(6) <E> 0(3) symmetry. Then we present issues related to the 
photoexcitation helicity amplitudes and the partial decay widths of nucleon resonances. 



A. Chiral constituent quark model 



As in Ref. 15j we start from an effective chiral Lagrangian 161 ]. 



(1) 
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where vector (V 1 ) and axial {A^) currents read, 

V » = l ( ^t + ^ d ^) f A m = i_ ( ^t _ ( 2 ) 

with £ = exp (i(f> m / f m ) and / m the meson decay constant. ?/> and m are the quark and 
meson fields, respectively. 

There are four components for the photoproduction of pseudoscalar mesons based on the 
QCD Lagrangian, 

(NflH^iNjlHjNi) 



+ M T , (3) 



where Ni(Nf) is the initial (final) state of the nucleon, and u(u m ) represents the energy of 
incoming (outgoing) photons (mesons). The first term in Eq. ([3]) is a seagull term. It is 
generated by the gauge transformation of the axial vector in the QCD Lagrangian. This 
term, being proportional to the electric charge of the outgoing mesons, does not contribute 
to the production of the ^-meson. The second and third terms correspond to the s- and 
it-channels, respectively. The last term is the t-channel contribution. 

In this paper we focus on the nucleon resonance contributions. Given that the tt-channel 
contributions are less sensitive to the details of resonances structure than those in the s— 
channel, it is then reasonable to treat the it-channel components as de gene rate 



For s-channel, the amplitudes are given by the following expression 
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2M N * kW 
M N * = — 2 N e -^<D N *, 4 

where yfs = W = E N + u 7 = E s + u m is the total centre-of-mass energy of the system, 
and On* is determined by the structure of each resonance. T(q) in Eq. (T4J) is the total width 
of the resonance, and a function of the final state momentum q. 
The transition amplitude for the n th harmonic-oscillator shell is 

O n = Ol + 3 n . (5) 

The first (second) term represents the process in which the incoming photon and outgoing 
meson, are absorbed and emitted by the same (different) quark. 
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In the present work, we use the standard multipole expansion of the CGLN ampli- 



tudes 



32j, and obtain the partial wave amplitudes of resonance hi,2i±i- Then, the transition 



amplitude takes the following form: 



O 



A'* 



*/i{±0" • e + f2i±& • qcr ■ (k x e) + if 3 i±a • kq • e + ifu±a ■ qe • q. 



(6) 



Expressing the CGLN amplitudes in their usual formulation 



33 



341 ]. leads to the Hebb- 



Walker amplitudes in terms of photoexcitation helicity amplitudes, 

A l± = +fA» 
Bi± = ±/ 1 

where 



\ N* 

1/2) 



1(1 + 2 



A N* 



./ 



1 



f MjVjgjV /c] 1/2 



2M N » 



-AT/2 = fo 



2M N * 



(2J + 1)2tt l M 2 m ' UJ s-M%* +iM N *T(q)^ 1/2 J ° s - M%. + iM N *T(q) 
with A™ 2 the N* — > r^A" decay amplitude, appearing in the partial decay width, 



1 [qljW ljm 2 

(2J+l) 7 rM JV J Al/2/Gm "' 



V 



raAf > 



(7) 
(8) 

(9) 

(10) 



where C* N represents the Clebsch-Gordan coefficients related to the isospin coupling in the 
outgoing channel. 



In Ref. 



15| . the partial decay amplitudes are used to separate the contribution of the 



state with the same orbital angular momentum L. In fact, with the helicity amplitudes of 
photon transition and meson decay we can directly obtain the CGLN amplitudes for each 
resonances in terms of Legendre polynomials derivatives: 



fu± — fo[+Ay 2 



f2l± — fo[+Ai/2 — 

hi± — =t/o- 



1 + 1/2 + 1/2 N * 
\\l + l/2±3/2 A3 ^ l±1 ' 



1 + 1/2 ±3/2 N . 

^/ + 1/2T 1/2 3/2j * 

n A N* 
Z/1 3/2 



4Z± 



+fo 



(I - 1/2 ± 1/2)(Z + 3/2 ± 1/2) 

OAN* 
Z/1 3/2 



P" 



P" 



(11) 



/ (Z-l/2±l/2)(Z + 3/2± 1/2) 

All fiS are proportional to the meson decay amplitudes. So they can be used to separate 
the contributions from the state with the same A" and L as presented in Ref. 



1 ■'-•- ] 1 A as well as the 
decay amplitudes, are related to the matrix elements of the electromagnetic interaction 



In our approach, the photoexcitation helicity amplitudes Ay* 2 and A^ 2 , 



Hamiltonian 



A x 
A™ 



>?l{N*;J\\H e \N;±\-l), 
(N;^u\H m \N*;Ju). 



(12) 
(13) 



B. Configuration Mixing 



The amplitudes in Sec JII Al are derived under the SU(6) ® 0(3) symmetry. However, for 
physical states that symmetry is broken. An example is the violation of the Moorhouse 
rule [35|. In Ref. 17], a set of parameters CV* were hence introduced to take into account 
the breaking of that symmetry, via following substitution: 



O n . ^C n *O n *. (14) 

In Refs. [nj, [ji]], those parameters were allowed to vary around their SU(6) <8> 0(3) values 
(I Cat* | = or 1). In this work, instead of using those adjustable parameters, we introduce the 
breakdown of that symmetry through the configuration mixings of baryons wave functions. 
To achieve such an improvement, we must choose a potential model. The popular used 



ones are one-gluon-exchange (OG 



model 



391 ] . As shown in Refs. 



40 



) model 



36 



37 



and Goldstone boson exchange 



4l|, these two models give similar mixing angles for the 



negative parity resonances and the relevant observables. Here, we adopt the OGE model 
which has been successfully used to study the helicity amplitudes and decay widths [13j of 
resonances. 

In OGE model, the Hamiltonian of system can be written as 



36 



37J,l38[, 



3 3 2 3 3 

# = E™* + E^2+ E \ Kr l + E u( rij )+H hyp , 

■i i<j=l «<i=i 



i=l 



(15) 



where the irti is the " constituent" effective masse of quark i and r ^ = — tj the sepa- 
ration between two quarks. The confinement potential has two components; one written 
as a harmonic oscillator potential (^Krfj, with K the spring constant), and an unspecified 
anharmonicity Ufaj), treated as a perturbation. 
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The hyperfine part interaction is the sum of contact and tensor terms, 

Hny P = f£ £ { y S„ ■ Sj S^) + -^( 3S "y - S, • Si ) } . (16) 



Here, Sj is the spin of quark i, and a s a normalization factor, treated as free parameter [38 ]. 

The hyperfine interaction generates the configuration mixings among the ground-state 
N 2 S S ([56,0+]) and other configurations, e.g. N 2 S' S ([56', 0+]), N 2 S M ([70,0+]), and 
N A D M ([70,2+]). Here, the notation is X 2 ^ 1 ^, where X = N, A, S the total 

quark spin, L = S, P, D... the total orbital angular momentum, and it = S, M or A the 
permutational symmetry (symmetric, mixed symmetry, or antisymmetric, respectively) of 
spatial wave function. 

The first two terms in Eq. ffl5l) can be rewritten as two harmonic oscillators within 
the Jocabi coordinate. Its solution is the well known SU(6) (g> 0(3) wave functions. The 
breakdown of the symmetry arises from the additional terms. Given that the configuration 



mixing is mainly produced by the spin- and flavor-dependent parts of Hami 



tonian 



401 ]. here 



371 ] . where three 



we use a simple method to deal with the confinement terms in Refs. [19 . 
constants E , Q, and A are introduced. 

In order to illustrate the modifications of the scattering amplitudes due to the SU(G) ® 
0(3) symmetry breakdown, we give in the following the explicit derivations in the case of 



the Sii(1535) resonance . In lines with Ref. [la], we express the amplitudes As n in terms 



of the product of the photoexcitation and meson-decay transition amplitudes, 

^ Sll oc< N\H m \S n >< S n \H e \N >, (17) 

where H m and H e are the meson and photon transition operators, respectively. The wave 
function can be written within the SU(6) ® 0(3) symmetry for n < 2 shells as X 2S+1 L n J p 
and configuration mixings, with J p the state's total angular momentum and parity, 

|Sii(1535)) = -sin^|iV 4 P M |-)+cos^|iV 2 PMr). ( 18 ) 
\Nucleon) = c x \N 2 S s \ + ) + c 2 \N 2 S' s ± + ) + c 3 \N 4 D M ± + ) + 

c,\N 2 S A 4 + } + c 5 \N 2 P A ^ + } } (19) 

where 9s and Cj can be determined by the OGE model. If we set c\ = 1 and £2,3,4,5 = 
(so 63 = 0), then, the SU(6) ® 0(3) symmetry is restored. The improvement compared 
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to Ref. [18|] is that here we not only take into account the mixing in the intermediate Sn 
resonance but also in the initial- and final-state nucleon. Moreover, for other resonances, we 
also include directly the configuration mixing of wave functions via OGE model, so that we 
do not need to introduce the free parameters C^* (Eq. (|14p ). 

The electromagnetic transition amplitudes then take the following form: 

< S n \H e \N > = Cl < S n \H e \N 2 Ssl + ) + c 2 < S n \H e \N 2 ^± + > + c 3 < S lx \H e \N A D M \ + ) 
+ c 4 < S u \H e \N 2 S M \ + ) + c 5 < S 11 \H e \N 2 P A ± + ) 

= Cl cos9 < N 2 P M \~\H e \N 2 S s \ + ) + ... (20) 



Here, the term < N A P M \ \H e \N 2 Ss\ + ) vanishes because of the Moorhouse rule 35]. In 



Ref. [IS], the mixing angles are introduced also to give a nonzero value for contributions 
from the _Di3(1700) resonance, but the nucleon wave function includes only the n = part, 
that is, C\ = 1, 02,3,4,5 = 0. Moreover, the contribution of the D 15 (1675) (\N 4 D M ^ + ) state) 
is zero, if we consider only the wave function up to n = 2. Then, in Ref. [18j], for this latter 
resonance a term identical to the contribution to the rj photoproduction on neutron target 
was added by hands. In this work, the nucleon wave function with n = 2 produces naturally 
a non-zero contribution with the same form as for neutron target under the 577(6) <8> 0(3) 
symmetry. 

Analogously, for meson decay amplitudes we get, 

<N\H m \S u > = c l (cose s -nsme s )<N 2 S s ^\H m \N 2 P M f} + ... (21) 
and the ratio 

<iV|/7 m |iV( 4 PM)i- > 

<iV|/7 m |iV( 2 PM) 1 - >' 1 1 

2 

is a constant determined by the SU(6) <8> 0(3) symmetry. 
Then, Eq. ( |T7|) reads, 

A Su = C Sll < N 2 S s \ + \H m \N 2 P M \- >< N 2 P M \~\H e \N 2 S s \ + > +..., (23) 

where, 

C Sll = c?(cos 2 9 S - K sin 9 S cos 6 S ) + ... (24) 

So, if we remove all n = 2 parts from the wave function of the nucleon, as in Ref. 171 ] . then 
the factor Cs u is a constant. However after other contributions are included, it becomes 
dependent on the momenta k and q. In this work we keep this dependence. 
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III. RESULTS AND DISCUSSION 



With the formalism presented in SecJTTl we investigate the process 7J9 — > rjp. A chiral 
constituent quark model was proven [3] to be an appropriate approach to that end. That 
work embodied one free parameter per nucleon resonance, in order to take into account 
the breaking of the 577(6) <8> 0(3) symmetry. In the present work, this latter phenomenon 
is treated via configuration mixing, reducing the number of adjustable parameters. As in 
Refs. 18), we introduce resonances in n < 2 shells, to study the rj photoproduction in the 



centre-of-mass energy W < 2 GeV. 



A. Fitting procedure 

Using the CERN MINUIT code, we have fitted simultaneously the following data sets: 



Differential cross-section: Data base includes 1220 data points, for 1.49 < W < 
1.99 GeV, coming from the following labs: MAMI [l|, CLAS [2|, ELSA 
and GRAAL |5j. Only statistical uncertainties are used. 



B, LNS J 



Polarized beam asymmetry: Polarized beam asymmetries (184 data points), for 
1.49 < W < 1.92 GeV, from GRAAL Q and ELSA Q. Only statistical uncertainties 
are used. 

Spectrum of known resonances: For spectrum of known resonances, we use as 
input their PDG values [42| for masses and widths, with the uncertainties reported 
there plus an additional theoretical uncertainty of 15 MeV, as in Ref. [l4], in order to 
avoid overemphasis of the resonances with small errors. The data base contains all 12 
known nucleon resonances as in PDG, with M < 2 GeV, namely, 

n=l: 5n(1535), 5n(1650), L>i 3 (1520), D u (1700), and L> 15 (1675); 

n=2: Pn(1440), P u (1710), Pi 3 (1720), P u (1900), F 15 (1680), F 15 (2000), and 
Fi 7 (1990). 

Besides the above isospin-1/2 resonances, we fitted also the mass of A(1232) resonance. 
However, spin-3/2 resonances do not intervene in the rj photoproduction. 
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• Additional resonance: Resonances with masses above M ~ 2 GeV, treated as 
degenerate, are simulated by a single resonance, for which are left as adjustable pa- 
rameters the mass, the width, and the symmetry breaking coefficient. 

The adjustable parameters, listed in Table [H are as follows: r\ nucleon coupling {(Jt/Nn) > 
mass of the non-strange quarks {m q ), harmonic oscillator strength (a), QCD coupling con- 
stant (a s ), confinement constants (Eq, Q, and A), three parameters M, T, and related 
to the degenerate treatment of resonances with masses above ~ 2 GeV, and the strength of 
the Pi3(1720) resonance. We will come back to this latter parameter. 

The spectrum of the known resonances put constraints on six of the adjustable parame- 
ters. Five of them (m q , a, a s , Q, and A) are determined through an interplay between the 
mass spectrum of the resonances and the photoproduction data via the configurations mix- 
ings parameters q (Eq. [201) . The sixth one, Eq, is determined by the mass of nucleon. The 
coupling constant g V NN is determined by photoproduction data. The parameter Cp 13 (i 72 o) is 
the strength of the Pi3(1720) resonance, that we had to leave as a free parameter in order 
to avoid its too large contribution resulting from direct calculation. This latter parameter, 
as well as those defining the higher mass resonance (HM N*) are determined by the pho- 



42) 



toproduction data. Notice that in fitting the photoproduction data, we use the PDG 
values for masses and widths of resonances. 

The complete set of adjustable parameters mentioned above, leads to our model A (see 
3 rd column in Table [J for which the reduce d y 2 turns out to be large (12.37). 



In recent years, several authors 



0, Q. 



22 



23 



24 



25 



27 



28 



29 



30, 



3 II ] have put 



forward need for new resonances in interpreting various observables, with extracted masses 
roughly between 1.73 and 2.1 GeV. We have hence, investigated possible contributions from 
three of them: Su, Di 3 , and Di 5 . For each of those new resonances we introduce then 
three additional adjustable parameters per resonance: mass (M), width (r), and symmetry 
breaking coefficient (Cn*)- Fitting the same data base, we obtained a second model, called 
model B, for which the adjustable parameters are reported in the last column of Table [B 
The reduced \ 2 is very significantly improved going down from 12.37 to 2.31. In the rest of 
this Section, we concentrate on the model B. 

Extracted values within OGE model come out close to those used by Isgur and Cap- 
stick 3, 0: E = 1150 MeV, Q » 440 MeV, A « 440 MeV. For three other parameters, 



10 



TABLE I: 
MeV. 



Adjustable parameters and their extracted values, 



with m q , a, E , S7, A, M, and T in 



Parameter Model A Model B 





Clnr> AT AT 


0.391 


0.449 






277 


304 




(X 


288 


285 






1.581 


1.977 




l!jc\ 


1135 


1138 




n 


450 


442 




A 


460 


460 






0.382 


0.399 


HM N*\ 










M 


1979 


2129 




r 


124 


80 




( j AT* 


-0.85 


-0.70 


New Sii: 










M 




1717 




r 




217 




Cj AT* 




0.59 


New D13'. 










M 




1943 




r 




139 




Cat* 




-0.19 


New D 15 : 










M 




2090 




r 




328 




O/v* 




2.89 


Xd.o.f 




12.37 


2.31 



Isgur and Capstick introduce 5 = (4a s a)/(3\/2irml), for which they get 300 MeV. Model 
B gives 5 « 262 MeV. 



For the three new resonances, we follow the method in Ref. 17| . as discussed in Sec JIIB| 
via Eq. (fT4j) . The extracted Wigner mass and width, as well as the strength for those 
resonances are given in Table [B 



23 
26 



or the new Sn, the Wigner mass and width are consistent with the values in Refs. 



18 



31 



but the mass is lower, by about 100 to 200 MeV, than findings by other authors 



22, 
2J, 



271 . l28l . |43j. The most natural explanation would be that it is the first Sn state in 



the n = 3 shell, however its low mass could indicate a multiquark component, such as, a 



quasi-bound kaon-hyperon [22] or a pentaquark state [441 ] . For the -D 13 (1850), the variation 
of x 2 is small. Interestingly, we find large effect from a -D15 state around 2090 GeV with a 
Wigner width of 330 MeV. It is very similar to the A[f2070)-D 15 reported in Refs. {3], 9]. It 
can be explained as the first 7J 15 state in n = 3 shell [3j. 

The results of baryon spectrum extracted from the present work are reported in Tables [III 
and IIHI Table [TT] is devoted to the known resonances. Our results are in good agreement 



with those obtained by Isgur and Karl 



33, 



and except for the Sii(1535) and _Di 3 (1520), 



fall in the ranges estimated by PDG 42j. The additional "missing" resonances generated by 



the OGE model, are shown in 
reported by Isgur and Karl 37|, 



Die IIHI The extracted masses are compatible with those 



38]. 



TABLE II: Extracted masses for known resonances. For each resonance, results of the present 



work (M OGE ) are give n in the first line, predictions from Isgur and Karl for negative-parity [38| 
and positive-parity I37I excited baryons in the second line, and PDG values [2] in the third line. 





#11(1535) 5n(1650) 


Pn(1440) Pn(1710) Pi 3 (1720) Pi 3 (1900) 


M OGE 


1473 1620 


1428 1723 


1718 


1854 


Refs. [37, 38] 


1490 1655 


1405 1705 


1710 


1870 


M PDG 


1535 ± 10 1655±^ 


14401^ 1710 ± 30 


17201^0 


1900 


D 13 (1520) D 13 (1700) Z?i 5 (1675) Fi 5 (1680) Fi 5 (2000) Fi 7 (1990) 


M OGE 


1511 1699 


1632 1723 


2008 


1945 


Refs. [37, 38] 


1535 1745 


1670 1715 


2025 


1955 


M PDG 


1520 ± 5 1700 ±50 


1675 ± 5 1685 ± 5 


2000 


1990 



In Table HV] we examine the sensitivity of our model to its ingredients by switching off one 
resonance at a time and noting the x 2 , without further minimizations. As expected, the most 
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TABLE III: Predicted masses for "missing" negative parity excited baryon by the present work 
(M OGE ) and by Isgur and Karl [37] . 





Pn 


Pu 


Pl3 


P13 


Pl3 


*15 


M OGE 


1899 


2051 


1942 


1965 


2047 


1943 


Ref. [37] 


1890 


2055 


1955 


1980 


2060 


1955 



important role is played by the S , n(1535), and the effects of Sn(1650) and _Di3(1520) turn 
out to be very significant. Within the known resonances, the other two ones contributing 
largely enough are -F 15 (1680) and Pi3(1720). In addition to those five known resonances, a 
new Sn appears to be strongly needed by the data, while the smaller effect of a new D i5 is 
found non-negligible. Finally, higher mass resonance (M ^ 2 GeV) and a new -D13 do not 
bring in significant effects. 

TABLE IV: The x 2 s shown are the values after turning off the corresponding (known) resonance 
contribution within the model B, for which \ 2 = 2.31. 



Removed N* 


5n(1535) 


5n(1650) 


P n (1440) 


Pn(1710) 


Pis(1720) 


Pi 3 (1900) 


x 2 


162 


11.9 


2.29 


2.39 


4.15 


2.35 


Removed N* 


^13(1520) 


£i3(1700) 


£>i 5 (1675) 


Pi 5 (1680) 


Pi 5 (2000) 


Fit (1990) 


x 2 


9.83 


2.29 


2.24 


4.82 


2.33 


2.31 


Removed N* 


HM N* 


New Sn 


New D13 


New D15 






x 2 


2.50 


12.69 


2.63 


3.88 







Our model B is built upon resonances given in Table HVl In Table |V] we investigate 
possible contributions from the missing resonances (Table HLT1) . Here, we add them one by 
one to the model B, without further minimizations. As reported in Table IV] none of them 
play a noticable role in the reaction mechanism. Please notice that for those resonances we 
use the masses that we have determined. We have checked the changes of the x 2 by varying 
those masses by ±100 MeV. Moreover, given that there is no unique information available 
on their widths, we have let them vary between 100 MeV and 1 GeV. The effects of those 
procedures on the reported x 2 s in Table IVl come out to be less than 10%. 

After having discussed above the s-channel contribution, we end this Section with a few 
comments. In our models, non-resonant components include nucleon pole term, and u- 
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TABLE V: The x 2 s shown are the values after adding the corresponding (missing) resonance 

contribution within model B, for which x 2 = 2.31. 

Added N* Pn(1899) Pn(2051) Pi 3 (1942) Pi 3 (1965) Pi 3 (2047) Pi 5 (1943) 
v 2 2.31 2.31 2.26 2.31 2.32 2.28 



channel contributions, treated as degenerate to the harmonic oscillator shell n. t-channel 
contributions due to the p- and ^-exchanges [45] , found [46j to be negligible, are not include 
in the present work. Our finding about the effect of higher mass resonances being very small, 
supports the neglect of the t-channel, due to the duality hypothesis (see e.g. Refs. [rsl. I^). 

Finally, the target asymmetry (T) data 48j are not included in our data base. Actu- 
ally, those 50 data points bear too large uncertainties to put significant constraints on the 



parameters 



46|. 



B. Differential cross section and Beam asymmetry 

In Figures HJ El and [HI we report our results for angular distributions of differential cross 
sections, excitation functions, and polarized beam asymmetries (£), respectively. Results 
for the models A and B are shown on all three Figures. The first striking point is that 
model A compares satisfactorily with data up to W < 1.65 GeV, but shows very serious 
shortcomings above, especially in the range W ~ 1.7 GeV to 1.8 GeV. Model B reproduces 
the differential cross section and polarization data well enough, though some discrepancies 
appear at the highest energies and most forward angles (W ^ 1.85 and 9 ^5 50°). 

In Fig. [I} we concentrate on the role played by the three most relevant known resonances 
discussed in Sec JIII Al (see Table HVl) . namely, by removing one resonance at a time, within 
the model B. The 5*11(1535) is by far the most dominant resonance at lower energies and 
has sizeable effect up to W « 1.8 GeV, while the S , n(1650) shows significant contributions 
only at intermediate energies. The _Di 3 (1520) has less significant contribution, but its role 
is crucial in reproducing the correct shape of the differential cross section, especially at 
intermediate energies. 

The importance of the other two known resonances, leading to a significant increase of 
X 2 when switched off (see Table HV|) . are illustrated in the left panel of Fig. [2j While, the 
Pi 3 (1720) affects extreme angles around W ~ 1.8 GeV, the Pi 5 (1680) is visible only at 
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FIG. 1: Differential cross section for the process 7 p — > rj p. The curves are for models A 
(dash-dot-dotted) and B (full). The other curves are obtained within model B by switching off 
one resonances at a time: Sn(1535) (dashed), 5*11(1650) (dotted) and L>i3(1520) (dash-dotted). 
The data are from CLAS (squares) [2], ELSA (circles) [3j, Mainz (diamonds) [1], and GRAAL 

(stars) 



forward angle. 

The right panel of Fig. [2] is devoted to the roles played by the three new resonances. 
As mentioned above, the main shortcoming of the model A appears around W ~ 1.7 - 1.8 
GeV. This undesirable feature is cured in the model B, due mainly to the new Sn, the mass 
of which turns out to be M = 1.717 GeV. Fig. [2] illustrates the increase of x 2 (Table HV) 
when that resonance is switched off in the model B. Smaller contributions from the new 
D15 appear in the forward hemisphere, while the new -D13 has no significant manifestation. 

Polarized beam asymmetry results are reported in Fig. [31 As shown in the left panel 
of that figure, although the model B gives a better account of the data than the model 
A, the contrast is less important compared to the differential cross-section observable. 
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W (GeV) W (GeV) 



FIG. 2: Differential cross section as a function of W at three angles. The dash-dot-dotted and 
full curves correspond to the models A and B. All other curves are obtained within the model 
B by turning off one known resonance or adding a missing one. In the left panel switching off 
Pi 3 (1720) (dash-dotted), Fi B (1680) (dashed); adding Pi 3 (1942) (dotted). In the right panel: 
switching off £11 (1730] (dash-dotted), £>i 3 (1850) (dotted), L>i 5 (1950) (dashed). The data are 
from CLAS (squares) Q|, Mainz (diamonds) [1], LNS (uptriangles) jjj. 

The Sii(1535) continues playing a primordial role, while the effect of 5n(1650) tends to 
be marginal. This is also the case (middle panel) for the known Pi 3 (1720) and missing 
Pi 3 (1942). The established importance of the Di 3 (1520) and F 15 (1680) (in left and middle 
panels, respectively) within this observable appear clearly. 

In the right panel of Fig. [31 we examine the case of three new resonances. The new Su 
gives sizeable contributions in the energy range corresponding roughly to its mass. In con- 
trast to the differential cross-section, the new _D 13 appears to be significant in the backward 
hemisphere. Finally, switching off the new D15 improves the agreement with the data at 
most backward angles shown, while for the cross-section we get an opposite behavior. This 
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FIG. 3: Polarized beam asymmetry as a function of W . The curves in the left panel are as in 
Fig. [TJ and those in middle and right pannels as in Fig. [2j The data are from ELSA (full circles) (| 
and GRAAL (stars) [5]. 

isolated contradiction reflects the relative weight of data for the two observables (roughly 
6 times more differential cross-section data than polarization asymmetry, with comparable 
accuracies) . 

This Section, devoted to the observables of the the process 7 p — > rj p, in the energy range 
W ^2 GeV, leads to the conclusion that within our approach, the reaction mechanism is 
dominated by five known and two new nucleon resonances. 



C. Helicity amplitudes and partial decay width 

As discussed in Sec. IV (Eqs. (20), (21), and (25)), our approach allows calculating the 
helicity amplitudes and the partial decay width N* — > 77 A^ within a given model without 
further adjustable parameters. 
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TABLE VI: Helicity amplitudes and decay widths for resonances, with r^y G = Ttot ■ # r »?7V in PDG 
42]. 



Resonances 


A/2 


. p nri 

A l/2 


A 3/2 


^3/2 




(^)V r ™ 


5n(1535) 


72 


90 ± 30 






7.05 


(+)8.87±i;i? 


5n(1650) 


60 


53 ± 16 






-2.20 


. in Q/i 

1 Qv ay4 
i - yo -1.57 


#u(1440) 


37 


-65 ± 4 










#il(1710) 


27 


9 ± 22 






1.30 


9 4Q+ 1 - 75 


#11 


3 








-1.64 




#11 


-2 








-0.76 




Pi 3 (1720) 


194 


18 ± 30 


-72 


-19 ± 20 


2.07 


1 1 ClA 

z,oo -0.71 


P 13 (1900) 


33 




1 




-0.87 


^ ^ _. | 9 1 1 
o or-hZ.ll 
O.O0_ 2 20 


#13 


32 




-2 




1.80 




As 


14 




2 




0.05 




#13 


-4 




4 




-0.73 




£>i 3 (1520) 


-20 


-24 ± 9 


144 


166 ± 5 


0.30 


^ — - _i_n D7 

o.5i±8:gg 


#13(1700) 


-6 


-18 ± 13 


2 


-2 ± 24 


-0.57 


0.00lJ;S 


#15(1675) 


-6 


19 ± 8 


-9 


15 ± 9 


-1.74 


u - uu -o.oo 


#15(1680) 


14 


-15 ± 6 


142 


133 ± 12 


0.44 


o-oo±S:oo 


#15 


-12 




5 




0.78 




#15(2000) 


-1 




13 




-0.38 




Fi 7 (1990) 


6 


1 


8 


4 


-1.25 


o-ooi^ 7 



In Table I VI I we report on our results within the model B, for all n =1 and 2 shell 
resonances generated by the quark model and complemented with the OGE model. In that 
Table, 2 nd and 4 th columns embody our results for the helicity amplitudes. Those amplitudes 



are in lines with results from other similar approaches (see Tables I and II in Ref. 19(). 



Comparing our results for the dominant known resonances of the model B with values 



reported in PDG 42J (3 rd and 5 th columns in Table IVll) leads to following remarks: i) Au 2 
amplitudes for S , n(1535) and Sii(1650), as well as Ay 2 and A 3 / 2 for # 13 (1520) and A 3 / 2 
for # 15 (1680) are in good agreement with the PDG values. For this latter resonances the 
A x j2 has the right magnitude, but opposite sign with respect to the PDG value. However, 

18 



for that resonance A 3 / 2 being much larger than A 1( / 2 , the effect of this latter amplitude is 
not significant enough in computing the observables. The amplitudes for Pi 3 (1720) deviate 



significant 
II in Ref. 



y 



rom their PDG values, as it is the case in other relevant approaches (see Table 
19]). Those large values produced by our model forced us to leave the symmetry 
breaking coefficient for P 13 (1720) as a free parameter (Table [I]) in order to suppress its 
otherwise too large contribution. As much as other known resonances are concerned we 
get results compatible with the PDG values for Pi 3 (1700) and Pi7(1990), and to a less 
extent for P 15 (1675). For P n (1440) our result deviates significantly from the PDG value. 
Once again, our result confirms the general trend observed in other works (see Table II in 
Ref. [19]), which very likely reflects the still unknown structure of that resonance. Finally, 
we put forward predictions also for the missing resonances, for which we find rather small 
amplitudes, explaining the negligible roles played by them in our model. 

The 6 th and 7 th columns in Table IVII show our results and PDG values, respectively, 
for the partial decay widths of resonances decay in the rjN channel, where a is the sign 
for it N — > 7] N as in Ref. 13]. Notice that the sign (a) in the PDG is known only for 
Sii(1535). Except for the two star resonance Pi3(1900), the theoretical results are close to 
the PDG values. 

It is worthwhile noticing that all dominated resonances in our model B have large helicity 
amplitudes, while some of them turn out to have rather small decay widths to the rjN 
channel. This result indicates that in looking for appropriate reactions to search for missing 
resonances it is not enough to have rather sizeable decay width, but one needs to put forward 
predictions for the observables. 



IV. SUMMARY AND CONCLUSIONS 



A formalism bringing together a chiral constituent quark approach and one-gluon- 
exchange model was presented and used to derive photoexcitation helicity amplitudes and 
partial decay width of the nucleon resonances. 

Our approach gives a reasonable account of the measured observables for the process 
7 p — > 7] p from threshold up to W ~ 2 GeV. Among the twelve nucleon resonances in that 
energy range, compiled by PDG, five of them are found to play crucial roles in the reaction 
mechanism, namely, Sn(1535), Sn(1650), Pi 3 (1720), P>i 3 (1520), and Pi 5 (1680). However, 
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those known resonances led to our model A, which does not allow an acceptable description 
of the data. Five extra resonances generated by the formalism, and known as missing 
resonances, turn out to show no significant contributions to the process under investigation. 
However, two new resonances reported in the literature, Su and Di 5 , are found relevant to 
that process; the most important effect comes from the Su resonance. We extracted the 
mass and width of those resonances: S u [1.730 GeV, 217 MeV], and D 15 [2.090 GeV, 328 
MeV]. Our model B, embodying those latter resonances, describes successfully the data. 

The helicity amplitudes and decay widths are calculated with the same parameters. Our 
results are compatible with other findings and come out close to the PDG values in most 
cases. 

To go further, we are pursing our investigations in two directions, 

• In the present work the s-channel resonances with masses above 2 GeV were treated as 
degenerate, given that the transition amplitudes, translated into the standard CGLN 
amplitudes were restricted to harmonic oscillator shells n < 2. recently, we have 
extended our formalism and derived explicitly the amplitudes also for n= 3 to 6 shells. 
Model search, including all known one to four star resonances in PDG, for W ~ 2.6 



GeV is in progress [46|] 



Our constituent quark approach applied to the 7 p — > K + A channel 



49j , showed that 



the intermediate meson-baryon states, treated within a coup 
have significant effects on the photoproduction observables 



ed channel formalism 



50j, 



3JJ. A more sophisticated 
coupling-channel treatment [5lJ has been developed and is being applied to the rj 
photoproduction reaction. Results of that work will be reported elsewhere. 
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APPENDIX A: MIXING COEFFICIENTS OF THE WAVE FUNCTIONS 



In Table IVH[ we present the mixing coefficients of the wave functions. In Ref. 



36, 



371 ]. Isgur and Karl have given their explicit values for positive parity and negative parity 
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resonances respectively. But in Ref. 36| the mixing between n = and n = 2 shells is not 



considered. Such mixings for the ground state are given in Ref. [13] without the contribution 
of 2 Pa- The parameters in that reference are determined only by the mass spectrum. Here 
we give our results by fitting both the mass spectrum and the rj photoproduction observables. 
In calculation we follow the conventions in Ref. {l3| . 



The mixing coefficients reported here lead to mixin. 



in agreement with results from other authors 



ang 



52|, |53j, |54j, 



55]. 



es, Q s = -31.7° and D = 6.4° 
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TABLE VII: Mixing coefficients of the wave functions. 



state 


wave function( 25 ' +1 L 7r ) 


Sll 


2 Pm 

- 1 lVl 


4 Pm 








N(1535) 


-0 851 


526 








N(1650) 




u.ou -L 








Pll 




a Dm 






2 Sm 


N(938) 


0.941 


-0.043 


-0.002 


-0.260 


-0.211 


N(1440) 


0.268 


0.000 


0.000 


0.964 


0.006 


N(1710) 


0.175 


-0.343 


-0.071 


-0.054 


0.919 




-0.103 


-0.839 


-0.424 


0.031 


-0.324 


N(2100) 


n rn? 


D 491 


U.c/UO 


n m n 

U.U-LU 


U.UOU 


P13 


'-'b 


2 LW 


ism 


2 Pa 


u M 


N(1720) 


0.858 


-0.483 


0.023 


-0.003 


-0.176 


N(1900) 


n ^14 




-U.OUJ 


\J.\JcfO 


u.ooy 




-0.185 


-0.482 


0.606 


-0.333 


0.505 




S59 


686 


496 


-0 S87 


-0 065 




-0.059 


-0.096 


-0.502 


-0.854 


-0.073 


D13 


2 p 

1V1 


4p 








N(1520) 


-0.994 


-0.111 








N(1700) 


-0.111 


994 








D15 


4 Pm 

J lVl 










N(1675) 


1.000 










F15 


2 D S 




4 D M 






N(1680) 


0.883 


-0.469 


0.001 








-0.457 


-0.860 


-0.225 






N(2000) 


-0.107 


-0.198 


0.974 






F17 


A D M 










N(1990) 


1.000 
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